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We study the fluctuations of the amplitude (i.e. the Higgs-Anderson) mode in a superconducting 
system of coupled Dirac particles proposed as a model for possible surface or interface supercon¬ 
ductivity in rhombohedral graphite. We show that the absence of Fermi energy and vanishing of 
the excitation gap of the collective amplitude mode in the model leads to a large fluctuation con¬ 
tribution to thermodynamic quantities such as the heat capacity. As a consequence, the mean-field 
theory becomes inaccurate indicating that the interactions lead to a strongly correlated state. We 
also present a microscopic derivation of the Ginzburg-Landau theory corresponding to this model. 


I. INTRODUCTION 

A topological flat band spectrum can emerge as a sur¬ 
face or interface state of topological semimetalJi^. Such a 
spectrum has beenproposed to occur in rhombohedrally 
stacked graphite^lEl, at dislocation interfaces in bernal 
stacked graphit^, edges of zig-zag graphen^ and at the 
interfaces in topological crystalline insulator^!. Such a 
state has a singular density of states which can, in the 
presence of an attractive pairing between the particles, 
lead to superconductivity with an unusually high critical 
temperaturdS*2l. This mechanism has been suggested as 
an explanation to the unusually high critical temperature 
Tc in some IV-VI semiconductor heterostructureJi^®^ 
and in various graphite based materialJ^^^. In the Hub¬ 
bard model, related work has also been done to study 
superconductivity in the presence of a flat bancP^ and 
with Fermi energy close to Van Hove singularitieJi^. 

We study here the effect of fluctuations on the su¬ 
perconducting properties of flat band superconductors. 
We use the particular model proposed in Refs. |3] and |3] 
and also numerically analyzed in Ref. [TH| for rhombohe¬ 
drally stacked graphite. The model consists of N indi¬ 
vidual graphene layers with s-wave pairing between the 
Dirac electrons coupled to form a stack of rhombohedral 
graphite. The flat band is formed on the surfaces of the 
stack where the superconductivity also appears (see Fig. 
[^. A peculiarity of this model is the closing of the gap in 
the fermionic excitation spectrum at the limit of a large 
number of layers even with an isotropic s-wave mean- 
field order parameter as can be seen from the surface 
state spectrum (see Fig. [2^1 

Ei = {i-pyp%snAi+e^), ( 1 ) 

where = Ui\p/pfb\^, 7i is the interlayer coupling 
constant, pfb = hlvp is the width of the flat band 
and Aq is the mean-field order parameter. (We use ev¬ 
erywhere h = kg = 1.) This expression is valid for 
\p/pfb\ < 1 — 1/fV. Enclosing of the gap at the flat band 
edge is due to the fact that there the surface states pene¬ 
trate into the bulk and surface superconductivity cannot 
create a gap for the bulk excitations. 



FIG. 1. Schematic figure of the model which consists of N 
parts with Dirac spectrum H = vpa ■ p coupled via coupling 
strength 71 . This results in an effective low-energy theory de¬ 
scribed by the action (10 1 . The snperconductivity is localized 
at the surfaces where the Gooper pair flelcP^is given by Ai,]v. 


Even without making detailed calculations, simple ar¬ 
guments as for why fluctuations are expected to matter 
in this model can be made. Typically the magnitude 
of fluctuations is characterized by the Ginzburg number 
Gi ^ 1 — T*/Tc, which is related to the temperature 
T* above which the fluctuations dominate the thermo¬ 
dynamical properties. In a 2D BCS superconductor, the 
Ginzburg number is given by Gi ^ T^/Ep, where Ep is 
the Fermi energy. For the flat band, Ep = 0, leaving no 
other relevant energy scale available to be compared with 
Tc which suggest^ that Gi ^ I. 

We first consider a Ginzburg-Landau theory for this 
model close to Tc. It is of the form expected from sym¬ 
metry considerations with coupling coefficients that we 
derive from the microscopic model. Based on this theory, 
we find that the contribution to the heat capacity from 
amplitude fluctuations is larger than the mean-field heat 
capacity jump at the superconducting transition for a 
wide range of temperatures, resulting in a large Ginzburg 
number of Gi = 2/5 in accord with the dimensional ar¬ 
guments above. Then we consider the same microscopic 
model, but without making the mean-field approxima¬ 
tion. Rather, we employ functional integral calculus to 
derive the partition function for the system that preserves 
fluctuations around the mean field up to the Gaussian 
approximation. We show that the correction from the 
amplitude mode to the mean-field heat capacity is large 
even far below the mean-field critical temperature Tc. Fi- 
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FIG. 2. Fermionic excitation spectrum of the surface state 
for superconducting (blue line) and normal (black line) state. 
For finite A^, the snperconducting state has an excitation gap 
Ep^min which closes in the limit N ^ oo. 


nally, we numerically show that close to T = 0 the ampli¬ 
tude mode and the mean field contributions to the heat 
capacity are approximately equal to each other. This 
finding can be analytically understood by assuming that 
at low temperatures the heat capacity is determined by a 
free boson contribution corresponding to the amplitude 
mode dispersion and a fermionic contribution due to Bo- 
goliubov quasiparticles. 


II. GINZBURG-LANDAU PICTURE 


A simple picture of the fluctuation contribution can 
be obtained from Ginzburg-Landau theory close to T^. 
There, a superconductor is described in terms of the 
Cooper pair held A whose absolute value as well as the 
gradient are small so that an expansion with respect to 
them can be made. In the usual approach, only the 
value of the wave function that minimizes the Ginzburg- 
Landau free energy is considered, assuming that all the 
pairs are condensed to this state. The study of huctua- 
tions in superconductors involves considering also Cooper 
pairs that are on some other states. In the held integral 
formalism we use below, this translates into taking an 
integral over diherent conhgurations of the Cooper pair 
helds and weighting different conhgurations by a factor 
that is given by the Ginzburg-Landau action when cal¬ 
culating observables. 

From symmetry considerations, or from microscopic 
calculations (see below), we get for a general supercon¬ 
ductor consisting of two separate symmetric parts (in our 
case, the top and the bottom layers of the graphite stack), 


a Ginzburg-Landau action of the form 
Sgl = [ (fx 


(oilAip-|-Q!2|VAip-|- 

iGfl.Af} 


—7(AiAAr + AatAi) 


( 2 ) 


Here 7 is the Josephson coupling coefficient between the 
two parts characterized by amplitudes Ai and Ajv. N 
in the index is the number of layers and it is assumed to 
be large. The partition function, from which thermody¬ 
namic quantities can be calculated, is then given by 


2 = n (3) 


A microscopic calculation from the general action for our 
model, given in the next section, yields the following val¬ 
ues for the Ginzburg-Landau coefficients: 


py _ Pfb T-Te _ _ 1 _ 

“ 487I-T2 Ta ’ ^2 — 32ttT^ ’ 

f> _ Pfb ry _ 'Iip'fb 

^ ■ 


(4) 


From these coefficients one can extract a coherence length 
~ Q! 2 /|ai| = 3Tc/(2p|.^|Tc — T|), which is very short 
(^ ~ wf/ti vp/A) for T sufficiently below Tc indicat¬ 
ing the strong coupling nature of the superconductivity. 

To consider fluctuations for T < we expand around 
the mean-field value of Aq = — (oi — 7 )/( 2 / 3 ) as Ai{x) = 
(Aq -I- The fluctuation part of the action 

can then be written in terms of four fluctuation modes: 


51 = 

5 2 = 

5 3 = 

54 = 
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(5) 


so that Sgl = Smf + Si + S 2 + S 3 + S 4 , where the 
mean-held part is given by Smf = —{oti — 7)^/(2/3). 
The four fluctuation modes can be identihed as the to¬ 
tal phase (Nambu-Goldstone) mode (^i), the total am¬ 
plitude (Higgs-Anderson) mode {S 2 ), the relative phase 
(Leggett) mode (S 3 ), and the relative amplitude mode 
(£'4). In this new diagonal basis, the four independent 
helds are dehned as 4>± = (j)i± (j)N and i5± = Si zL 6 n- In 
the limit A —>■ 00 , the two phase and the two amplitude 
modes become identical and the action separates into two 
identical standard huctuating Ginzburg-Landau actions 
with coefficients given by Q. The standard Ginzburg- 
Landau theory calculations can then be applied^ and we 
hnd that the total amplitude mode contribution to the 
heat capacity below become^ 


sa 


amplitude 


^Pfb'^c 
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( 6 ) 
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Here A is the surface area of the sample. In the lowest 
order in Tc — T, the phase mode does not contribute to 
the heat capacity. Setting (5C'ampiitude(2^*) equal to the 
heat capacity jump at the transition, 


AC — Appg 


5 

Gtt 


J {T-T,Y ^ 


(7) 


and solving for temperature, T*, yields the Ginzburg 
number 



This is of the order of unity as expected from purely di¬ 
mensional arguments. While at the reduced temperature 
indicated by the Ginzburg number, the Ginzburg-Landau 
approach is not, strictly speaking, valid, this nevertheless 
gives us an estimate of the size of the fluctuations. 

To understand the large fluctuations and compare 
them to conventional superconductors, we note that 
in both cases the fluctuation contribution and the 
heat capacity jump are proportional to the same 
Ginzburg-Landau coefficients, but the microscopic val¬ 
ues of the coefficients differ greatly. In a conven¬ 
tional two-dimensional BGS superconductoil^ ai = 
N{0){T - T,)/Tl a 2 = 7a^)N{0)vl/{32 tt^T^) and 
/3 = 7C(3)fV(0)/(I67r^r^), where ({x) is the Riemann 
zeta function, A^(0) is the density of states at the Fermi 
level and uf is the Fermi velocity. There is one addi¬ 
tional parameter compared to the case of the flat band 
model. This leads to different scales for the heat ca¬ 
pacity jump and the fluctuation heat capacity in two- 
dimensional BGS superconductors 


deduce that the action has the form 


s = Y,^{p) 


_ p' %UJ' 


'p _ 

^N,p—p' 


'i'b') 
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p,p' ^ p 


( 10 ) 

Here Cp = 1i{\p\/pfb)^, Wp = (2n-h l)7rr/(l -pVpfb) 
is the Matsubara frequency with a momentum depen¬ 
dent factor which comes from using the ansatz for sur¬ 
face state wave functions given in Ref. [Sj and g is the 
superconducting coupling strength with dimensions of 
(energy) • (volume). The volume V = Ad, where d ^ 
1 /pfb is the spatial extent of the surface state wavefunc- 
tions. Here, p = (n, p) is the three-momentum which 
includes the two-dimensional momentum p in the plane 
of the layers and the Matsubara index n. The 4x4 ma¬ 
trix structure of the action (in addition to the momentum 
space degrees of freedom) results from top-bottom layer 
G particle-hole degrees of freedom. The Gooper pairs are 
described by matrices 


A, 


,p-p> 


f _ 0 ^i,p—p' j 

\\y_p 0 J 


( 11 ) 


and can now have a non-zero momentum, i.e., p — p' ^ 0. 
Integrating over the fermionic modes in (101 results in 


Z= PAPA exp 


-^5](iAi,pp-f iA^,pn 

. ^ p 


+trln 


( 12 ) 
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Tc Tc 
Ef \T-Tc\ 


AC 


(9) 


resulting in a much smaller Ginzburg number Gi ^ 
Tc/Ep compared to the flat band superconductors. 


III. MICROSCOPIC CALCULATION AND 
ARBITRARY TEMPERATURES 


At low temperatures, the Ginzburg-Landau approach 
is invalid and we need a microscopic theory to account 
for the fluctuations. For the surface state, we can use 
the Bogoliubov - de Gennes equation derived in Ref. [3] to 


where the trace tr is taken over the three-momentum and 
matrix indices. 

From expression (12), the Ginzburg-Landau action, 
([^, with the coefficients can be obtained by ex¬ 
panding in small A^ p. The mean-field value for the 
Cooper pair field Aq(T) is found by minimizing the ac¬ 
tion and Tc is found by solving for the highest tem¬ 
perature below which the mean-field value is non-zero. 
In the limit N —>■ oo, this yields the mean-field value 
Ao(0) = pp’^pQ/ilddir) and Tc = Ao(0)/3 which we have 
used in expressing the Ginzburg-Landau coefficients in 
terms of the critical temperature. More generally, we can 
expand A.^ around this mean-field value, A^ = Ag -I- SAi. 
This separates the action into a mean-held part and a 
huctuation part, S = Smf + Ss, where the fluctuation 
part is given by 
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(13) 


and the mean-field part by 

Smf = -45]ln [cosh(£;p/2T)] + —IAq^. (14) 

p ^ 

In the latter we have performed the summation over the 
Matsubara frequencies and are the positive quasipar¬ 
ticle energies given by Eq. Q. The partition function is 
then given as an integral over the fluctuation modes as 

Z = [ VsAvsA^e-^^. (15) 


We have here written the action in a convenient ma¬ 
trix form with the vector fluctuation fields given by 
6 Aq = {6Ai^q,d~Ai-q,6AN^q,S~AN-q)^- The matrix ele¬ 
ments in the fluctuation action are given in terms of three 
polarization operators. Expressions for them are 


p 

p 

'^. = E 


UJpLJp—q 

i^p + ^0 + Q)i^p-q + ^0 + ^p-q) 

^0 

i^p + ^0 + Q)(.^p-q + ^0 + ^p-q) 
i^p + ^0 + Q){^l-q + ^0 + ip-q) 


(16) 


As for the Ginzburg-Landau action, the four eigenmodes 
of the fluctuation action correspond to relative and total 
phase fluctuations and relative and total amplitude fluc¬ 
tuations. Also, again, in the limit of a large number of 
layers, only the total amplitude and total phase modes 
survive and the third polarization operator, Xq, vanishes. 

Erom the Ginzburg-Landau calculation, we know that 
the dominating contribution to the thermodynamics in 
the large N limit comes from the amplitude fluctuation 
part. Taking into account only this mode, the fluctua¬ 
tion contribution to the free energy, 6 F = —TlnSZ « 
^Eanipiitude: becomes 


SE 


amplitude 


2TEln 1 


Wq-Vq\ 

V/gT ) 


(17) 


We calculate ^Eampiitude as well as the heat capacity, 
iJC'ampiitude = -T’^f.^Fampiitude, numerically. The re¬ 
sult for the free energy is shown in Eig. Eor zero 
temperature, the fluctuation contribution goes to zero, 
but grows larger than the mean-field contribution as the 
temperature increases. Above the critical temperature, 
i^Eampiitude Saturates into a finite value. 

Eor the heat capacity, the numerical result is shown in 
Fig.S At low temperatures we find that the fluctuation 



T/T„ 


FIG. 3. Free energy contribution from the amplitnde fluctu¬ 
ations and from the mean-field calculation. 



FIG. 4. Specific heat contribution from the amplitude fluc¬ 
tuations (green curve) and from the mean-field calculation 
(blue curve). For comparison the fluctuation contribution for 
a conventional BGS superconductor with Tc = 10~‘^£'_f and 
equal heat capacity jump at the transition AC is also shown 
(red curve). Inset shows the low-temperature behaviour of 
the three curves. 


contribution to the heat capacity becomes equal to the 
mean field heat capacity Cmf, both given by 


C{T « E) = 


Apfb^ T 
6 An’ 


(18) 


For Cmf Eq. (18) straightforwardly follows from 
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Eq. (14). On the other hand, the fact that (JCampiitude is 
also given by Eq. (18) can be analytically understood by 


assuming that it comes from a free boson contribution 
corresponding to the amplitude mode dispersion, which 
can be determined from the corresponding eigenvalue of 
the matrix in the fluctuation action (13). Namely, this 
way we find that for T = 0 and N ^ oo the amplitude 
mode dispersion is given by 


Ep,e.nipi{T = 0) = AoP^/pI-b, 


(19) 


and the corresponding heat capacity is indeed given by 
Eq. (18). For higher temperatures, the fluctuation contri¬ 


bution completely dominates over the mean field contri¬ 
bution. At Tc, the heat capacity diverges as ^ 1/\T — Tc\ 
as predicted by the Ginzburg-Landau theory. The tem¬ 
perature at which the fluctuation heat capacity equals 
the mean-field heat capacity jump at the transition is 
approximately T /Tc ~ 0.4 which is on the same order as 
expected from the Ginzburg number. For comparison, we 
also show the Ginzburg-Landau fluctuation contribution 
expected for a conventional BGS superconductor [Eq. (|^] 
with Tc/Ep ~ 10“^ and other parameters chosen so that 
the mean-field heat capacity jump at the transition AC 
would be the same as the one we get for the flat band 
case. 


IV. DISCUSSION 


The large contribution from the amplitude mode even 
at low temperatures can be understood from the form 
of the amplitude mode dispersion at T = 0. By inves¬ 
tigating the eigenvalues of the matrix in the fluctuation 
action, Eq. (13), one can show that the mass gap of the 


mode is not 2Ao as in conventional superconductors, but 
instead it is 2 Ep,mim where Ep^rnin is the minimum of the 
fermionic spectrum. Since the spectrum of the fermionic 
excitations is given by Eq. (1^, which does not have a 
gap in the limit of large N ^ee Fig. [^, also the gap 
in the amplitude mode closes and the contribution from 
the fermionic mode and the bosonic amplitude mode can 
be of the same order of magnitude at low temperatures. 
Therefore, finite N corrections to our result would likely 
diminish the fluctuation contribution, improving the va¬ 
lidity range of the mean-field theory. 


Besides rhombohedral graphite, there have been also 
other suggestions for correlated states at the surfaces of 
topological semimetals. We expect that their mean-field 
spectrum is also gapless in the bulk limit, and proba¬ 
bly exhibits very similar fluctuation contributions as that 
described in this paper. Moreover, flat bands can also 
emerge at the surface in superconductors with a nodal 
order parameter structure possibly resulting in correlated 
states^^’^. Similarly as in the present situation, the 
fluctuation contribution to the heat capacity may also 
dominate the mean field contribution in these correlated 
states, so that the mean field theory becomes inaccurate 
and the interactions lead to a strongly correlated state. 

We consider here only the effect of the amplitude mode 
to thermodynamics since we assume that it dominates 
over the phase mode contribution. At reasonably high 
temperatures it follows from general Ginzburg-Landau 
theory that the amplitude mode contribution dominates 
because of its singularity at T = Tc. At lower temper¬ 
atures, the phase mode contribution to thermodynamic 
properties in a charged system is a lso in fluenced by the 
coupling to the electromagnetic fielcPSEH. In the calcula¬ 
tion of this coupling it should be taken into account that 
although the order parameter appears only close to the 
surface the supercurrent flows also in the bullP^. 

Lastly, we want to emphasize that we do not consider 
the coupling of the fluctuations to external fields. The 
large contribution from the amplitude fluctuations to the 
heat capacity might or might not appear in other observ¬ 
ables depending on whether they couple to this degree of 
freedom or not. For example, it is known that the elec¬ 
tromagnetic field does not couple to the amplitude mode 
at T = 0 directlji^ and thus such a large fluctuation 
effect is not necessarily present in the current. 
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